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Abstract 



The density matrix of a qudit may be reconstructed with optimal efficiency if 
the expectation values of a specific set of observables are known. In dimension six, 

■ the required observables only exist if it is possible to identify six mutually unbi- 
ased complex (6 x 6) Hadamard matrices. Prescribing a first Hadamard matrix, 
we construct all others mutually unbiased to it, using algebraic computations per- 
formed by a computer program. We repeat this calculation many times, sampling 
all known complex Hadamard matrices, and we never find more than tzuo that are 

^ , mutually unbiased. This result adds considerable support to the conjecture that no 

■ seven mutually unbiased bases exist in dimension six. 

^3 ; 

1 Introduction 

Suppose you want to reconstruct the density matrix p of a qudit, a quantum system 
i/-^ I with d orthogonal states. To apply the most efficient reconstruction method you should 

O ' measure a set of observables associated with d mutually unbiased complex Hadamard 

matrices of size {dx d). A complex Hadamard matrix if is a unitary matrix having entries 
of modulus only; two such matrices are said to be mutually unbiased (MU) if their 
^ ' product is another Hadamard matrix, 

H^H' = H" , (1) 

^ I where denotes the adjoint of the matrix H. The columns of d MU complex Hada- 

■ mard matrices and of the identity define a complete set of 1) MU orthonormal bases 
in the space suitable for optimal state reconstruction [HEI. Complete MU bases are 
often characterized directly by the scalar products between their elements. 



5jj, if 6 = 6' , 
73 if^^^'' 



where b,h' = 0, 1, . . . , d; they are also used in quantum cryptography [3] and play an 
important role in the solution of the Mean King's problem [4J. 

Here is the catch: as of today, complete sets of MU bases have been constructed 
only in spaces of prime or prime power dimension. If the dimension is a composite 
number, d = 6, 10, 12, . . ., the existence of a complete set of MU bases in C"' has neither 
been proved nor disproved (see [5] for a review). In other words, it is unknown even 
for a qubit-qutrit system whether there exists a set of observables which would realize 
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optimal state reconstruction. Interestingly, constructing complete sets of MU bases is 
equivalent to finding an orthogonal decomposition of the Lie algebra s/rf(C), which 
poses a long-standing open problem whenever d is not a prime power [6J. 

Let us summarize what is know about the (non-) existence of MU bases in compos- 
ite dimensions. There are a few analytic results: 

• it is possible to construct three MU bases in without reference to the value of d 
IZi; hence, three MU bases do exist for any composite dimension d; 

• there are at least {p^ + 1) MU bases if p, the smallest factor of the prime decom- 
position of d, occurs k > 1 times 18J; 

• more than {p^ + 1) MU bases are known exist for specific values of d; for example, 
if d = 2^ X 13^, a total of 6 (= + 2) MU bases have been identified f9l. 

Attempts to generalize number-theoretic formula used in the construction of complete 
MU bases from prime-power dimensions to composite dimensions fail [lOj. Further- 
more, searches for MU bases in dimension six by numerical means have been unsuc- 
cessful: 

• no evidence for the existence oifour MU bases in has been found |TT|; 

• strong numerical evidence against the existence of various MU constellations (cor- 
responding to subsets of four MU bases) has been obtained, making the existence 
of a complete set highly unlikely fVT\. 

Some rigorous results have been obtained by restricting the search to MU bases of a 
specific form: 

• selecting a first Hadamard matrix and then searching for MU vectors with com- 
ponents given by suitable roots of unity leads to no more than two MU complex 
Hadamard matrices, or three MU bases in ffT3|: 

• Grassl |iZJ has shown that only finitely many vectors exist which are MU with 
respect to the identity and a given complex Hadamard matrix related to the 
Heisenberg-Weyl group in C^. Again, no more than two MU Hadamard matrices 
emerge, giving rise to at most three MU bases; it is thus impossible to base the 
construction of a complete set on the Heisenberg-Weyl group. 

The strategy of this paper will be to generalize Grassl's approach by removing the 
restriction that the second MU basis be related to the Heisenberg-Weyl group. Instead, 
we will consider many different choices for the second MU basis, thoroughly sampling 
the set of currently known complex Hadamard matrices. We will find that none of the 
matrices studied can be used to construct a complete set of MU bases. Taken together, 
these negative instances provide further strong support for the conjecture that no seven 
MU bases exist in dimension six. 

Let us now present the outline of our argument. In Sec. |2l we briefly describe the 
set of known complex Hadamard matrices in dimension six. Then, we explain in Sec. 
|3] how to construct all vectors that are MU with respect to both the standard basis of 

and a second basis, defined by an arbitrary fixed Hadamard matrix. We illustrate 
the algorithm for d = 3 only to rediscover the known complete set of four MU bases. 
Then, whilst rederiving Grassl's result for d = 6, we will explain the subtle interplay 
between algebraic and numerical calculations in this approach. Sec. S] presents our 
findings which we obtain by applying the algorithm to nearly 6000 Hadamard matrices 
of dimension six. Conclusions are drawn in the final section. 
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2 Complex Hadamard matrices in dimension six 



Traditionally, a Hadamard matrix H in dimension d is understood to have elements ±1 
only and to satisfy the condition H^^H = dl, where / is the identity. In the context of 
MU bases, it is customary to call H a Hadamard matrix if it is unitary and its matrix 
elements are of the form 

\H,,\ = ^, i,j = 0,1,... ,d-l. (3) 

The d vectors formed by the columns of such a matrix provide an orthonormal basis 
of C^. Each of these vectors is mutually unbiased with respect to the standard basis, 
naturally associated with the identity matrix /. It is convenient to identify a Hadamard 
matrix with the MU basis formed by its columns. 

Two Hadamard matrices are equivalent to each other, H' w H, if one can be obtained 
from the other by permutations of its columns and its rows, and by the multiplication of 
its columns and rows with individual phase factors. Explicitly, the equivalence relation 
reads 

H' = M1HM2 , (4) 

where Mi and M2 are monomial matrices, i.e. they are unitary and have only one 
nonzero element in each row and column. Consequently, each Hadamard matrix is 
equivalent to a dephased Hadamard matrix, the first row and column of which have 
entries 1 / Vd only. 

All (complex) Hadamard matrices are known for dimensions d < 5 but there is 
no exhaustive classification for d = Q. It is useful to briefly describe the Hadamard 
matrices known to exist in dimension six since we will 'parametrize' the search for 
MU bases in terms of Hadamard matrices. We use the notation introduced in |[l4i| the 
authors of which maintain an online catalog of Hadamard matrices L15J . 

Each point in Fig. [H an updated version of a figure presented in |T3|, corresponds to 
one Hadamard matrix of dimension six, except for the interior of the upper circle where 
a point represents two Hadamard matrices (cf. below). There is one isolated point, rep- 
resenting the spectral matrix S given in [16J, also known as Tao's matrix [17J. Three sets 
of Hadamard matrices labeled by a single parameter are known: the Dip family D{x) in- 
troduced in HH, a family of symmetric matrices denoted by M{t) [19] and the family of 
all Hermitean Hadamard matrices B{9) [20 J. Two two-parameter families of Hadamard 
matrices are known to arise from discrete Fourier-type transformations F{xi,X2) in C^, 
and from their transpositions, -F^(xi, X2) l|2T| . The SzoUosi family X{a, h) is the only 
other known two-parameter set |i22| . Interestingly, the matrix X{Q, 0) can be shown to 
be equivalent to F(l/6, 0), and there is a second possibility to define a matrix at this 
point, giving rise to -'^^(0, 0) F^(l/6, 0) [23J. We have noticed that such a doubling 
actually occurs for all values of the parameters (a, h) leading to a set of Hadamard ma- 
trices X^{a, b) inequivalent to X{a, h). Hence, the interior of the upper circle in Fig. [T] 
represents two layers of Hadamard matrices which are glued together at the boundary 
of the circle. Topologically, the Szollosi family X{a, b) and the set X^ {a, b) thus com- 
bine to form the surface of a sphere. Appendix |A] lists the explicit forms of Hadamard 
matrices as well as the parameter ranges which have been reduced to their fundamental 
regions using the equivalence relation (jH). 

Fig. [1] also shows equivalences between Hadamard matrices simultaneously be- 
longing to different families. The circulant Hadamard matrix C ||24|, for example, em- 
beds into the Hermitean family which in turn is given by the boundary of the Szollosi 
families. Lining up some of the points where different families overlap suggests that 
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we arrange the Hadamard matrices in a symmetrical way. There are two reasons for 
this: for the Di^a, Hermitean and symmetric families, a reflection about the line passing 
through the points F(0, 0) and S maps -ff (x) to x), while the same reflection sends 
i?(x) to i7^(x) if the matrix H{x.) is a member of the Di^a, Hermitean or Fourier fami- 
lies; for the SzoUosi family, the reflection about the vertical axis must be supplemented 
by a change of layer in order to get from X{a, b) to X'^{a, h). 

Let us finally mention that the known families of Hadamard matrices come in two 
different types, affine and non-affine ones. The set -ff (x) is affine if it can be written in 
the form 

H{x) = H{0) o Exp[i?(x)] (5) 

for some matrix R; the open circle denotes the Hadamard (elementwise) product of 
two matrices, [A o B)ij = AijBij, and Exp[i?] represents the matrix R elementwise 
exponentiated: (Exp[R])ij = exp Rij. Both Fourier-type families and the Di^a matrices 
are affine (cf . Appendix |A]) while the symmetric, Hermitean and SzoUosi families are 
not. 



3 Constructing MU Vectors 

In this section, we make explicit the conditions on a vector to be MU with respect 
to the standard matrix and a fixed Hadamard matrix, i.e. tothe pair {/, H}. Then we 
outline an algorithm to construct all solutions of the resulting multivariate polynomial 
equations, allowing us to check how many additional MU Hadamard matrices do exist. 
We illustrate this approach by constructing a complete set of four MU bases in dimen- 
sion d = 3, and we reproduce Grassl's result for d = 6 in order to explain that this 
approach produces rigorous results in spite of inevitable numerical approximations. 



3.1 MU vectors and multivariate polynomial equations 

A vector \v) G is MU with respect to the standard basis (associated with the columns 
of the identity /) if each of its components has modulus 1 / Vd. Furthermore, \v) is MU 
with respect to a fixed Hadamard matrix i7 if I = l/d, where \h{k)) is the state 

associated with the k^^ column h{k) of H, k = 0, . . . ,d — 1. 

Let us express these conditions on |u) in terms of its components vj, written as 



Vdvj 



1 

Xj + iyj 



l,...,d-l, 



(6) 



where Xj,yj are 2{d — 1) real parameters. The overall phase of the state \v) is irrelevant 
which allows us to fix the phase of its first component. Then, the first set of constraints 
on the state \ v) reads 



2 , 2 



1 



J 



1, 



,d-l. 



and the second set is given by 



(7) 



Eh*{k)^ 

j=0 



d-l 



^Hl.{xj+iyj) 



1 

d 



0,...,d-2, 



(8) 



where the state \h{k)) has components hj{k) = Hjk, = 1, . . . , d — 1. The completeness 
relation of the orthonormal basis {\h{k)), k = 0, . . . ,d — I}, implies that if a state \v) is 
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MU with respect to {d — 1) of its members, it is also MU with respect to the remaining 
one. Therefore, it is not necessary to include k = d — lin Eqs. (H)). 

For each given Hadamard matrix H, Eqs. (O and (H)) represent 2{d—l) simultaneous 
coupled quadratic equations for 2{d — 1) real variables. Once we know all solutions of 
these equations, we know all vectors \v) MU with respect to the chosen pair of bases 
{/, H}. Analysing the set of solutions will reveal whether they form additional MU 
Hadamard matrices, or, equivalently, MU bases. 

If Eqs. ([Zl and (H) were linear, one could apply Gaussian elimination to bring them 
into 'triangular' form. The resulting equations would have the same solutions as the 
original ones but the solutions could be obtained easily by successively solving for the 
unknowns. 

The solutions of Eqs. (O and (H) can be found using Buchberger's algorithm II25I 
which generalizes Gaussian elimination to (nonlinear) multivariate polynomial equations. 
In this approach, a set of polynomials V = {p„(x), n = 1, . . . , N} is transformed into 
a different set of polynomials Q = {gm{'^),m = 1, . . . , M} (usually with M ^ N) 
such that the equations V = and Q = possess the same solutions; here "P = is 
short for Pn(x) = 0,n = 1,...,A^. Technically, one constructs a Grobner basis Q of the 
polynomials V which requires a choice of ordering the variables [25 J. The transformed 
equations Q = Q will be straightforward to solve due their 'triangular' form: one can 
find all possible values of a first unknown by solving for the zeros of a polynomial in a 
single variable; using each of these solutions will reduce one or more of the remaining 
equations to single-variable polynomials, allowing one to solve for a second unknown; 
etc. This process iteratively finds all solutions of ^ = and, therefore, all solutions of 
the original set of equations, V = 0. 

A Grobner basis exists for any set of polynomial equations with a finite number of 
variables. However, the number of steps required to construct a Grobner basis tends to 
be large even for polynomials of low degrees and a small number of unknowns. Thus, 
Buchberger's algorithm is most conveniently applied by means of algebraic software 
programs. We have used the implementation [27 \ of this algorithm suitable for the 
computational algebra system Maple [28J since we found it to be particularly fast for 
the system of equations under study. 

Let us now illustrate how to construct all vectors MU with respect to a pair {/, H} 
by solving the multivariate polynomial equations (|7| and ^ using Buchberger's al- 
gorithm. We will consider two cases in dimensions d = 3 and d = 6, respectively, 
which have been solved before but they are suitable to illustrate the construction and 
to discuss some of its subleties. 

3.2 Four MU bases in 

In dimension d = 3, four MU bases are known to exist. We will now show how to 
construct two MU Hadamard matrices H2 and H3 given a pair {/, H}. The resulting 
three MU Hadamard matrices plus the identity are equivalent to a complete set of four 
MU bases in C^. 

1. Choose a Hadamard In dimension three, all Hadamard matrices are known and 
there is only one choice for a dephased Hadamard matrix [21 J given by the Fourier 
matrix. 




(9) 
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where oj = exp(27ri/3) is a third root of unity. 



2. List the constraints We want to find all states \v) € C'^ which are MU with respect 
to the columns of the identity matrix / and the Fourier matrix F3. Using the four real 
parameters xi, X2, yi, and y2 introduced in the constraints (O and (H)) read explicitly 

l-xl-yl = 0, 

xi + X2 + X1X2 + yiy2 = 0, 
X1 + X2- \/3yi + \/3y2 + xiX2 - \/^xiy2 + V3yiX2 + j/ij/2 = . (10) 

The solutions of these four coupled quadratic equations in four real variables, V = 0, 
will tell us whether additional Hadamard matrices exist which are MU with respect to 
the Fourier matrix F3. 



3. Construct the solutions By running Buchberger's algorithm, we find the Grobner 
basis Q associated with the polynomials in Eqs. ([TO)) . Equating the resulting four poly- 
nomials gri(x), n = 1, . . . , 4, to zero, gives rise to the equations 

3y2-42/i = 0, 

l-X2-2yi = 0, 

1 + 2x1 + 4yiy2 - %i = 0, 

3 - 4y2 + 4yiy2 - 4yi = 0. (11) 

This set is 'triangular' in the sense that solutions can be found by iteratively deter- 
mining the roots of polynomials for single variables only. The first equation has three 
solutions, 

y2e{0,±\/3/2}; (12) 
next, the second equation implies that 

/ if y2 = , 
'' = [2 ify2 = ±V3/2; (13) 

etc. Altogether, there are six solutions, 

= -1, VS, VS), Sb = 2, -VS, 0) , 



Sr. 



i(2,-l,0,-V3), Sd = i(-l,-l,-\/3,-V3) 



Se = i(2,-l,0,\/3), s/ = i(-l,2,\/3,0), 

defining s = (xi, X2, yi, ^2)- 

Since the degrees of the polynomials Q in Eqs. (fTTj) does not exceed three, we are 
able to obtain analytic expressions for its solutions. This, however, is a fortunate coindi- 
cence due to the simplicity of the problem: in general, we will need to determine the 
roots of higher-order polynomials (cf . the example presented in Sec. I3.3|l which requires 
numerical methods. The resulting complications will be discussed in Sec. [ 
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4. List all MU vectors Upon substituting the solutions to sj into (|6]), one obtains 
six vectors 




which are MU with respect to the columns of both the matrices / and F3. No other 
vectors with this property exist, leaving us with Va, ■ ■ ■ ,Vf, as the only candidates for 
the columns of additional MU Hadamard matrices. 



5. Analyse the vectors The six vectors in ((14)) allow us to define an additional Hada- 
mard matrix only if any three of them are orthogonal; for a second Hadamard matrix 
the remaining three must be orthogonal among themselves and MU to the first three. 
Calculating the inner products between all pairs of the vectors Vatovj shows that they 
indeed fall into two groups with the required properties. Consequently, we have con- 
structed a complete set of four MU bases in C^, corresponding to the set {I, Fs,H2, H^} 
where the columns of the matrices H2 and H3 are given hy {va,Vb,Vc} and {vd,Ve,Vf}, 
respectively. 

We have also checked that the construction procedure works in dimensions d = 
2, 4, 5 and d = 7 where it correctly produces complete sets of d + I MU bases. The 
matrices F2, F3 and F5 are the only dephased Hadamard matrices in dimensions d = 2,3 
and d = 5, and there is only one way to construct complete MU bases from the vectors 
obtained. We have thus shown that the Heisenberg-Weyl construction of a complete 
set of MU bases is essentially unique in dimensions two, three and five. 



3.3 Three MU bases in 



In d = 6, the existence of seven MU bases is an open problem. We will search for all 
states \ v) which are MU with respect to the identity / and the six-dimensional equiva- 
lent of F3 given in ([9]), the dephased Fourier matrix 



/ 1 
1 
1 
1 
1 

V 1 



1 


1 


1 


1 


1 


UJ 










^2 




1 






LO^ 


1 




1 




LO^ 




1 




a;2 






Cj3 




UJ 



(15) 



with Lo = exp(7rz/3) now being the sixth root of unity. This problem has been studied 
in the context of biunimodular sequences [24J and in relation to MU bases [7J. It is im- 
possible to complement the pair {I,Fq} by more than one Hadamard matrix MU with 
respect to Fq. Thus, the construction method of MU bases in prime-power dimensions, 
which is based on the Heisenberg-Weyl group, has no equivalent in the composite di- 
mension d = 6. We will now reproduce this negative result. 

Having chosen the first Hadamard matrix to be Fq, we can write down the condi- 
tions which the components of a state \v) must satisfy, V = 0. After some algebraic 
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operations detailed in Appendix|Bl one obtains the equations 



xi + x^ + xiX2 + X2X'j, + X'iXi + XiX^ + yiy2 + y2y3, + ViyA + ym = 0, 
yi - + xiy2 - X2yi + X2y:i - x:iy2 + x^yi - XiV^. + x^yr^ - xr^y^ = 0, 

X3 + X1X4 + X2X5 + yiy4 + ^2^5 = 0, 
X2 + X4 + X1X3 + X1X5 + X2X4 + X3X5 + yiys + yiys + ^2^4 + 2/3^5 = 0, 
y2 - y4. + xiys - xiy5 + X2y4. - X3yi + x^yrj - X4y2 + x^yi - x^y-i = 0, (16) 

which must be supplemented by the five conditions (0 arising for d = 6. 

We need to find all solutions of these ten coupled equations V = which are 
quadratic in ten real variables. The Grobner basis G associated with the set V con- 
sists of 36 polynomials of considerably higher degrees. We reproduce only the first one 
of the new set of equations, G = 0, 

- 245025 2/5 + 4318758 y^^ - 28135161 ys^ + 89685000 2/5 - 158611892 y^^ 

+177275680 2/5 " - 150745472 7/5 + 104333824 - 43667456 2/5 (17) 
+23511042/5^^ + 48824322/5^^- 17039362/5^^ + 2621442/5^^ = 0, 

being of order 25 in the single variable 2/5. This equation admits 15 real solutions, 

2/5 G {0, ±1, ±^(1 + ±^(1 - V^), ±0.988940 . . . , ±0.622915 . . .} , (18) 

the last four of which we only find numerically. Due to the triangular structure result- 
ing from Buchberger's algorithm, there will be equations (at least one) containing only 
2/5 and one other single variable. For each value of 2/5 taken from ((T8|) , they reduce to 
single-variable polynomials the roots of which can be determined to desired numerical 
accuracy; etc. Keeping track of all possible branches we obtain 48 vectors that satisfy 
the Eqs. ((T6ll . 

Having determined the candidates for columns of MU Hadamard matrices, we cal- 
culate the inner products among all pairs of the 48 vectors. It turns out that there are 
16 different ways to group them into bases of C^. However, no two of these bases are 
MU with respect to each other. Consequently, it is possible to form at most 16 differ- 
ent triples of MU bases which include Fq, but it also follows that the Fourier matrix Fq 
(or any other unitarily equivalent element of the Heisenberg-Weyl group [7J) cannot be 
part of three MU Hadamard matrices giving rise to four MU bases. 

There are, however, many choices other than H = Fq for a dephased Hadamard 
matrix in dimension six. In Sec. IH we will repeat the calculations just presented for a 
large sample of currently known Hadamard matrices. Before doing so, we will discuss 
the fact that we are able to construct the desired vectors only approximately. In the 
following section we show that sufficiently high numerical accuracy allows us to draw 
rigorous conclusions about the properties of the exact vectors. 



3.4 The Impact of Numerical Approximations 

The previous section illustrated that the problem of finding MU vectors with respect 
to the identity I and a given Hadamard matrix H can be reduced to successively solv- 
ing for the roots of polynomials of a single variable. These roots, however, can only 
be found approximately. Does the approximation prevent us from drawing rigorous 
conclusions about the properties of the MU vectors we construct? We will argue now 
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that it remains possible to find upper bounds on the number of MU vectors with the 
desired properties. 

Consider the system of polynomials V = {p„(x), n = 1, . . . , 10} in the variables 
X gM^'^ resulting from some chosen Hadamard matrix H, and calculate a Grobner ba- 
sis, Q = {gm{^),Tn = 1, . . . , A/}. The roots of the equations "P = and Q = Q are 
identical by construction. Since Q = Q corresponds to a 'triangular' set, its roots can 
be found iteratively but, in general, no closed form will exist. The implementation of 
Buchberger's algorithm which we have chosen finds these roots with user-specified 
accuracy, relying on the theory presented in [29 J. 

Now suppose that Q = Q has two roots and s;,, to which we have found approxi- 
mations, sa and s^. The associated states \va) and \vb) will be close to the exact states, 
\va) and \vi,). It is convenient to think of the vectors \va) and \ vb) as lying inside a ball 
centered at the exact states, the radius of the ball being determined by the accuracy of 
the approximate roots. If the inner product of the exact states \va) and \vi)) has a non- 
zero modulus, 5 > 0, then they are not orthogonal. We can detect this by calculating 
the inner product of the approximate states \ va) and \vb)- If we ensure that the error of 
the approximate inner product is small compared to the number 5, a non-zero lower 
bound for the exact scalar product can be obtained, and it signals non-orthogonality of 
the exact states. A similar argument allows us to exclude that two states are MU with 
respect to each other. 

We determine the roots of ^ = to 20 significant digits which proves sufficient to 
put relevant limits on the properties of the vectors constructed in dimension six. The 
results presented in the next section thus provide strict limits on complementing pairs 
of bases {/, H} by MU vectors. 

4 Constructing MU Bases in Dimension Six 

We are now in a position to present the main results of this paper. We will consider one 
Hadamard matrix H a\.a time constructing all additional Hadamard matrices MU with 
respect to the chosen one. Picking matrices both systematically and randomly, we will 
find that not a single one is compatible with the existence of four MU bases. 

More specifically, we will calculate three quantities for each chosen Hadamard ma- 
trix H: the number which equals the number of vectors MU with the pair {/, H}; 
the number Nt providing an upper bound of the number of different triples of MU 
bases {I, H, H'}; and the value of Np, an upper bound of the number of constellations 
{5, 5, 1, 1} (or {/, H, \va), \vb)}), i-e. we report how many pairs of MU vectors can be 
found among those already MU with respect to / and H. 

4.1 Special Hadamard Matrices 

To begin, we consider Hadamard matrices which are either isolated or special in the 
sense that they belong to different families simultaneously. All these matrices are to 
be found on the symmetry axis of Fig. Q: the Fourier matrix Fq = F(0, 0) being 
invariant under transposition, the Di\a matrix Dq = D{0) which is both symmetric and 
Hermitean, the circulant matrix C, and the Spectral matrix S. 

The first row of Table [T] completes the findings of Sec. 13.31 obtained for the Fourier 
matrix Fq: there are Ny = 48 vectors MU with respect to both / and Fq, forming Np = 
144 MU pairs. The 48 vectors can be arranged in Nt = 16 different ways to form a 
second Hadamard matrix H' which is MU with respect to Fq. However, no two of 



9 



these 16 Hadamard matrices are MU between themselves, limiting the number of MU 
bases containing Fg to three. 

A similar analysis for the Dila matrix Dq reveals that there are 120 vectors MU to its 
columns and those of the identity, 60 of which form ten bases but none of these are MU 
with respect to each other. Whilst ten triples of MU bases exist, sets of four MU bases 
which include Dq do not exist since A'j, = 0. 

Interestingly, the components of the 120 vectors have phases cj) which take values in 
a small set only, 

(t)D = {0, TT, ±7r/12, . . . , ±ll7r/12, ±a} , (19) 

where tana = 2. This resullQ agrees with the one obtained by Bengtsson et al. lfT3| 
(note, however, that the descriptions given in the last two entries of the list in their Sec. 
7 must be swapped). What is more, our approach proves that these authors have been 
able to identify all vectors MU with the pair {/, Dq} by means of their ansatz for the 
form of MU vectors. In fact, the values of Nt and Np in Table [l] given for Dq are exact, 
not upper bounds, since the phases of the MU states \v) are known in closed form. 

The circulant matrix C permits only 38 MU vectors, none of which are MU with 
respect to each other, A'j, = 0, and none of their subsets forms a MU basis as A'^t = 0. 

The spectral matrix S is the only known isolated Hadamard matrix. We find 90 MU 
vectors but not a single sextuple of orthonormal ones among them while there are 1,115 
pairs of MU vectors. Thus, the pair {/, 5} cannot even be extended to a triple of MU 
bases. 

4.2 Affine Families 

Table |2] collects the properties of vectors MU with respect to the pair {I, H} where H is 
an affine Hadamard matrix, i.e. taken either from the one-parameter set discovered by 
Di^a or from the two-parameter Fourier or SzoUosi families. Again, we have sampled 
the relevant parameter spaces both systematically and randomly. 

The set of Di^a matrices D{x) depends on a single continuous parameter x, with 
\x\ < 1/8. We have sampled the interval in steps of size 1/144 making sure that the 
resulting grid of points include the 24th roots of unity which play an important role for 
Dq, so 

Fi) = {a/144 : o = ±1,±2...,±18}; (20) 

note that the matrix Dq has been left out. The number of vectors MU with the pair 
{/, D{x)} depends on the value of the parameter x: the Di^a matrices D{x) on the grid 
F D allow for either 72 or 120 MU vectors which can be grouped into into four additional 
Hadamard matrices. Since they are not MU between themselves, there are at most three 
MU bases containing any of these Di^a matrices. No MU pairs exist, Np = — except 
for two small intervals near x = ±1/8 where either 24 or 48 pairs can be found. 

The results obtained from randomly picking points in the fundamental interval are 
in line with the observations made for grid points. Fig. |2] shows N^, the number of 
vectors MU with respect to the pair {/, Z)(x)} for all 536 values of the parameter x 
which we have considered. The function N^{x) appears to be symmetric about x = 
and piecewise constant, dropping from 120 for small values of x to 72 at x ~ ±0.0177, 
and to 48 at the end points of the interval, x = ±1/8. We are not able to explain why 
Nv{x) = Ny{—x) should hold although we do know how to derive a similar property 

^Such a restricted set of phases also occurs for other members of the Di^a family. For example, all 48 
vectors MU with the pair {J, D(l/8)} have phases limited to the set (f>D U {±/3} where tan/3 — 3. 
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for the families of symmetric and Hermitean Hadamard matrices (cf. Sec. I4.3|l . The 
values for Np and can be found in Table HI 

The results for members of Fourier family F(x) are qualitatively similar. Picking 
values of X = (xi,X2) either randomly in the fundamental area or from the two- 
dimensional grid 

rp = {(a,6)/144:a = l,2, ...,24, 6 = 0, 1, . . . , 12, a > 26} , (21) 

invariably leads to 48 vectors being MU to the columns of the pair {/, F(x)}. There 
are eight different ways to form additional Hadamard matrices for each point consid- 
ered except for the matrix F(l/6, 0) with an upper bound of 70 triples. The number of 
pairs varies, Np G {0, 6, 12, 24}, without any recognizable systematic behaviour. It is 
important to realize that Grassl's result — the construction of complete sets of MU bases 
cannot be based on the Heisenberg-Weyl group in dimension d = 6 — also holds for the 
2,168 other Fourier matrices we have considered. 

The situation is similar when turning to the family of transposed Fourier matrices, 
F^(x). The number equals 48 throughout and a second Hadamard matrix can 
be formed in eight different ways, and only matrix F^(l/6, 0) allows for 70 different 
triples, eight being the norm. The number Np runs through the same values as in the 
Fourier family. 

4.3 Non-Affine Families 

The equations V = encoding MU vectors for the symmetric M{t), Hermitian B{9) 
and SzoUosi X{a,b) families turn out to be more challenging from a computational 
perspective: the program has, in general, not been able to construct the associated 
Grobner bases Q. The problem is not a fundamental one — the desired Grobner bases 
do exist but it appears that their construction requires more memory than the 16GB 
available to us. 

We suspect that the difficulties are due to the fact that, for non-affine matrices, the 
coefficients of the polynomials V = are no longer equal to fractions or simple roots 
of integers. When approximating the coefficients in question by fractions we obtain 
different sets of polynomials, V, and, interestingly, the program indeed succeeds in 
constructing the corresponding Grobner bases, Q, outputting (approximate) MU vec- 
tors \v). Being continuous functions of the coefficients, the approximate vectors will 
resemble the exact ones, \v) ~ \v). However, the number of MU vectors may change 
discontinuously if "P = is considered instead of P = 0, similar to the discontinuous 
change in the number N^ for the family D{x) near x ~ 0.0177, shown in Fig. |2l In other 
words, it could happen that we 'lose' some solutions due to a geometric instability as a 
consequence of modifying the defining polynomials. 

To determine the impact of such an approximation, we have studied how the num- 
ber Ny of MU vectors changes in a case for which we know rigorous bounds. We retain 
only five significant digits of the coefficients in the equations "P = associated with 
the family D{t) and solve for the approximate MU vectors. The inset of Fig. |2] shows 
that the plateaus of 120 and 72 MU vectors continue to be well-defined away from the 
discontinuity at a; ~ 0.0177 while the values of N^ fluctuate close to it. Assuming that a 
qualitatively similar behaviour will also occur for symmetric and Hermitean matrices, 
we now simplify the equations "P = associated with them. Retaining only five sig- 
nificant digits of the coefficients in these equations, we determine the number of MU 
vectors \v) and their inner products. 
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Fig. |3] shows that the family of symmetric Hadamard matrices M{t) comes with 
48 MU vectors \v) close to the point t = 0, while there are 120 near t = 1/4. These 
numbers are consistent with the rigorous bounds obtained in Sec. 14.21 if we recall that 
M(0) = M(l/2) F(0,0) and M(l/4) D{0) holds (cf. Fig. [H). Across the entire 
parameter range, the number of MU vectors is a piecewise constant function symmetric 
about X = 1/4, with distinct plateaus of 48, 52, 120 and possibly 96 MU vectors. We 
suspect that the other values of Ny near the discontinuities are spurious. An analysis 
of the scalar products among the approximate MU vectors shows that they cannot be 
arranged into a single additional basis. Table |3] lists the results obtained for both a 
regular grid 

Tm = {a/144 : a = 1, 2, . . . , 71; a / 36} ; (22) 

and 300 randomly selected points in the fundamental interval; the reason for leaving 
out a = 36 is the equivalence M(l/4) « Dq just mentioned. We are confident that a 
more rigorous approach will confirm the absence of triples of MU bases containing a 
single symmetric Hadamard matrix M{t). 

Further support for the reliability of the approximation follows from the overall 
symmetry of A'^j, in the fundamental interval. To show that must be an even fimction 
of the parameter t, let us make explicit the transformation of a symmetric matrix M{t) 
under complex conjugation, 

M{-t) = PM*{t)P , (23) 

where P is a permutation matrix. Now consider a vector \v) which is MU to the 
columns m{t) of the matrix M{t); then 

\{m{t)\v)\^ = \{m{t)\vT\' = \{m*{t)\v*)\^ = \{m{-t)P\v*)f , (24) 

and, therefore, P\v*) is MU with respect to each column of M{—t). 

Thus, the vectors MU to M{—t) are the complex conjugates of those MU to M{t) 
with entries permuted by P. Hence the number of vectors, Ny, is an even function of t 
and the vectors MU to M{t) have the same properties as those MU to M {—t). Although 
we did not pay any attention to the existence of this exact symmetries when introducing 
the approximations, the results obtained do respect it. 

The results obtained for Hermitean Hadamard matrices B{9), shown in Fig. |3l are 
similar to those of the symmetric family. The observed plateaus blend in with the rig- 
orous bounds found for Ny = 120 and A'^ = 38 due to the equivalences i?(l/2) w D(0) 
and 5(6*0) ~ C (cf. Tabled]). We consider the plateaus at 56, 58, 60, 72, 84 and 108 to be 
genuine while spurious values for Ny proliferate near their ends, where A'^j, is likely to 
vary discontinuously. Once more. Table |3] reveals that neither regularly spaced points 
on the grid 

Tb = {a/144 : a = 55, 56, . . . , 89; a / 72} ; (25) 

nor randomly chosen values of 9 define Hadamard matrices B{9) which would allow 
the construction of more than three MU bases. As for the symmetric family, we are 
able to explain the symmetry of Ny-. members of the Hermitean family satisfy the re- 
lation B{1 — 6) = B*{0) which implies that the number Ny of MU vectors (and their 
properties) will not change upon a reflection about the point 6 = 1/2. 

Finally, let us consider the Szollosi family, the non-affine two-parameter set of Hada- 
mard matrices X{a, b). Fig. |4]shows the values of Ny for randomly chosen parameters 
on two cuts through parameter space, namely along the line 

A = {(a,6) : arg(a + i6) =7r/6} (26) 
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which connects X(0, 0) « F{l/6, 0) to the circulant matrix C, and the randomly chosen 
line 

A' = {{a, b) : arg(a + ib) = 0.3510} (27) 

connecting X{0,0) to B{9'), a Hermitean Hadamard matrix on the boundary. The val- 
ues of Ny at the end points of the lines are, in both cases, consistent with results ob- 
tained above for F(l/6,0), C, and B{6'); broadly speaking, the number of solutions 
again represents a step function. However, the plateaus at 48, 52, 54, 56, 58 and 60 in 
Fig. m (b) show considerable overlap: the effect of approximating the coefficients in 
the relevant polynomials is even more pronounced for the SzoUosi family than for the 
other non-affine families. The results for the 300 randomly chosen parameter values 
sampling the two-dimensional parameter space resemble those of the symmetric and 
Hermitean families: we find 48 < A^t, < 120 throughout, and none of the selected 
Hadamard matrices allow for a triple of MU bases. Preliminary calculations show that 
the properties of the new family of transposed SzoUosi matrices X^(a, b) are similar to 
those of the set X{a, b). 

While not being exact, the results for the symmetric, Hermitean and SzoUosi fami- 
lies provide bounds on the number of MU bases which can be constructed from their 
members. Except in some special cases with higher symmetry, it is not possible to even 
find a third MU basis. We consider it unlikely that the approximation made would 
systematically suppress other MU vectors with properties invalidating this conclusion. 

5 Summary and Conclusions 

We have searched for MU bases related to pairs {/, H} where / is the unit matrix and 
H runs through a discrete subset of known (6 x 6) complex Hadamard matrices. Using 
Buchberger's algorithm, we have obtained upper bounds on the number of MU bases; 
the bounds are rigorous in many cases and approximate in others. Each of the 5,980 
calculations required between 4 and 16 GB of memory and, altogether, would have 
lasted approximately 29,000 hours on a single 2.2 GHz processor. 

Each point in Fig. |5] represents one of the Hadamard matrices H we have been 
investigating. We find that for any set of three MU bases, {/, H, H'}, at least one of 
the MU Hadamard matrices H and H' must either be a member of the Di|a or Fourier 
families. None of the remaining Hadamard matrices shown in Fig. |5]can be a member 
of a triple of MU bases. Further, any set of four (seven) MU bases in dimension six 
would require three (six) MU Hadamard matrices different from the ones shown in 
Fig. |5l This clearly conforms with the numerically obtained evidence that no four MU 
basis do exist [12] . 

There is one caveat that we must make regarding the results of the non-affine fami- 
lies. In general, the program was unable to construct the associated Grobner bases for 
the symmetric, Hermitian and SzoUosi families. For these Hadamard matrices, we are 
unable to guarantee that we have found all MU vectors, however, we consider it un- 
likely that the approximation made would systematically suppress the missing vectors. 

It has been suggested that the set of Hadamard matrices in depend on four 
parameters [13], a conjecture which recently gained some numerical support [30]. It 
remains difficult to draw general conclusions about the number of MU bases in dimen- 
sion d = 6. However, we would like to point out that the approach presented here 
is future-proof: it will work for any Hadamard matrix — including currently unknown 
ones. 
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In summary, we have shown that the construction of more than three MU bases in 
is not possible in a considerable number of instances. The results add significant 
weight to the conjecture that a complete set of seven MU bases does not exist in di- 
mension six. It becomes ever more likely that only prime-power dimensions allow for 
optimal state reconstruction. 
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Figures 



Spectral: S • 



Circulant: 

c^B(e„) 



Hermitean: 6(0) 



SzoUosi: X(a,b), 



D(l/8) 



F(l/6,0) 




D(-l/8) 



Dita: D(x) 



F'(l/6,0) 



Fourier: 



F(l/6,1/12) 



Fourier 
Transpose: 



F^(l/6,1/12) 



Figure 1: The set of known Hadamard matrices in dimension six consists of special 
Hadamard matrices F(0, 0) = Fq, D{0) = Dq,C, and S, located on the vertical symme- 
try axis; of the afftne families D{x), F(x), and -F^(x); and of the non-ajfine families M{t), 
B{9), X{a, b), and X^(a, b) (see Appendix |A] for definitions). Note that the sets X{a, b) 
and X'^{a, b) cover the interior of the upper circle twice. 
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-1/8 




X 



1/8 



Figure 2: The number of vectors | v) which are MU with respect to the columns of 
the identity / and Di^a matrices D{x); the parameter x assumes 72 parameter values 
X G r^), and 500 randomly chosen ones in the fundamental interval [—1/8, 1/8] of the 
parameter x. The inset illustrates the impact on Ny near the discontinuity x ~ 0.0177 if 
an ap-proximate set of equations is used (cf. Sec. I4.3|) . 
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(a) 



(b) 



120 



100 



A^„ 80 



60 



40 



Figure 3: The number Ny of vectors \v) which are MU with respect to the columns of 
the identity / and (a) symmetric Hadamard matrices M{t); the parameter t assumes 60 
parameter values t G Fm, and 300 randomly chosen ones in the fundamental interval 
[0, 1/2], and of (b) Hermitean matrices B{9); the parameter 6 assumes 34 parameter 
values 6 £ Fb, and 300 randomly chosen ones in the fundamental interval [9o,l — Oq]. 
The phase has been defined in equation ((35)) . 
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(a) 



56- 



• ••Ml 



54 



N 52 



50 



48 ••• • 



1/2 



\a+ib\ 



\a+ib\. 



(b) 62 




1 \a+ibl 



Figure 4: The number N.^ of vectors which are MU with respect to the columns of the 
identity / and Szollosi Hadamard matrices X{a, h) for 50 randomly chosen parameter 
values (a) on the line A connecting F(l/6, 0) to C, and (b) on the line A' connecting 
F(l/6, 0) to B{6'); in both figures, the maximum modulus \a + ib\ is defined by Eq. 
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Figure 5: The set of all Hadamard matrices H which have been considered (cf. Fig. [T] 
and Tables [nil): given {I, H}, only the Fourier matrices -F(x), the (transposed) Fourier 
matrices F^(x), and the Di\a matrices D{x) allow for the construction of a second MU 
Hadamard matrix and, thus, a third MU basis; in all other cases not even three MU 
bases exist. 
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Tables 



H 




Nt 


Np 


Fe 


48 


16 


144 


Do 


120 


10 





C 


38 








S 


90 





1115 



Table 1: The number of MU vectors and their properties for special Hadamard matrices: 
there are Ny vectors being MU with respect to the pair of matrices {/, H}; these vectors 
form Nt additional Hadamard matrices (i.e. there are A^^ different triples of MU bases) 
and Np is the number of pairs of vectors MU to themselves and to {/, H}. 
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#(x) 
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0/24/48 




random 


500 
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0/24/48 


F(x) 




168 


48 


8/70 


0/6/12/24 




random 


2,000 


48 


8 


0/6/12/24 


F^'(x) 


Tf 


168 


48 


8/70 


0/6/12/24 




random 


2,000 


48 


8 


0/6/12/24 



Table 2: The number of MU vectors and their properties for affine Hadamard matrices: 
the second column indicates which values have been chosen for the parameters x; the 
grids of points Tm and Tf are defined in Eqs. ((20]l and ((2T|) . respectively; the third 
column displays the number of Hadamard matrices considered in a sample; Ny, Nt 
and Np are defined as in Tableland vary as a function of the parameter values (cf. Sec. 
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#(x) 


N, 


Nt 


Np 


M{t) 


Tm 


70 


48-120 





0-18 




random 


300 


48-120 





0-18 


B{9) 


Tb 


34 


56-120 





0-14 




random 


300 


56-120 





0-14 


X{a,b) 


A 


50 


48/56 





0-10 




A' 


50 


48-60 





0-10 




random 


300 


48-120 





0-10 



Table 3: The number of MU vectors and their properties for non-affine Hadamard ma- 
trices: the grids Tm and Tb are defined in Eqs. (|22)) and (|25l) , respectively; see Eqs. ((26)) 
and ((27)) for the definition of the lines A and A'; other notation as in Table|2) preliminary 
results for the family X^(a, b) resemble those obtained for X{a,b). 
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A Known complex Hadamards matrices in dimension six 



This Appendix lists the currently known complex Hadamard matrices for easy refer- 
ence and to establish notation. For more details the reader is referred to |(T3| and to the 
online catalogue iTSf . 



A.l Special Hadamard matrices 

The Fourier matrix Fq has been introduced in Eq. ((15)) : it is contained in both the Fourier 
family -F(x) and the transposed Fourier family -F^(x) for x = 0, where Fq = F{0, 0) f« 
F'^iO, 0) holds (cf. Sec. lA^l) . 

The Difa matrix Dq is an example of a complex symmetric Hadamard matrix. 
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(28) 



embedded in a continuous one-parameter set of Hadamard matrices, the Di^a family 
(cf.lAl). 

Bjorck's circiilant matrix ||24| is defined by 
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(29) 



where 




(30) 



It was originally thought to be isolated but it is now known to be part of the family of 
Hermitian Hadamard matrices, C « B{9q) (cf. IA.3|l . 

The only known isolated Hadamard matrix is the spectral matrix, 
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(31) 



where w is a third root of unity, uj = e^'^'/^. It has been discovered by Moorhouse | [T6| 
and, independently, by Tao |fT7| . 
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A.2 Affine Families 



There are three affine families of Hadamard matrices, characterized by the property 
^ that they can be written as a non-trivial Hadamard product. The Dila family [18] is 
given by D{x) = Dq o Exp[27r iR{x)], \x\ < 1/8, with Dq from Eq. ((28)) and 
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(32) 



the componentwise exponential Exp[ ] of a matrix has been defined after Eq. 

The Fourier matrix Fq has been embedded in a similar way into a too-parameter 
set, namely the Fourier family -F(x) = Fq o Exp[27rii?(x)], where 



R{x) = R{xi,X2) 



/OO 000 0\ 

Xi X2 Xi X2 



Xi X2 Xi X2 



V Xi X2 Xi X2 / 



(33) 



the parameters (xi,X2) take values in a fundamental region given by a triangle with 
vertices (0,0), (1/6,0) and (1/6,1/12). 

Upon transposing the matrices F(x) one obtains a different two-parameter set of 
Hadamard matrices, called the transposed Fourier family F^(x). It has the same funda- 
mental region as the Fourier family. 

A.3 Non-Affine Families 

Non-affine Hadamard matrices cannot be written in the form ((5)). The Hermitean family 
pO| provides a one-parameter example of such a set. 
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(34) 



where y = e and t 



xyz, with 

l + 2y 



y(-l + 2y + y2) 



l + 2y + y^± ^2(1 + 2y + 2y^ + y^) 
1 + 2y - y2 

the free parameter 9 is restricted to vary within the fundamental interval [9q,1 — 9q\ 
and the number 9q is defined by the condition 
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Note that this is a smaller fundamental region than was previously known; the re- 
duction is due to equivalences that have become apparent since the discovery of the 
SzoUosi family (cf. below). 

Another non-affine one-parameter set of Hadamard matrices is given by the sym- 
metric family l[T9|. 
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(36) 



where x = e , and the complex numbers a, b, c, d,p, q are the unique solutions of the 
equations 

1 + x + d + a + b + c = 0, 
x"^ -2x -2a-2d-l = , 
1—x + b + c + p + q = 0, 

x^ + 2b + 2c+l = 0. (37) 



In addition, one needs the fact that given a row (ri, . . . , rg) of a Hadamard matrix, the 
last two elements are determined by E = (ri + r2 + + ri)/2, since 



(38) 



if S 7^ 0. The fundamental region is given hy t G [0, 1/2]. 
Finally, there is the non-affine SzoUosi family 1122 II 



X{a,b) = H{x,y,u,v) = 
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(39) 



0, respectively. 



The entries x, y and u, v are solutions to the equations fa = and /_c 
where 

fa{z) = z^-az'^ + a*z-l, (40) 
and a = a -\- ibis restricted to the region B defined by D{a) < and D{—a) < 0, with 

D{a) = \a\^ + 18|ap - 8Re[a^] - 27 . (41) 

It is possible to reduce D to a smaller fundamental region [23] since, firstly, the trans- 
formation a —a maps Hadamard matrices to equivalent ones and, second, Eq. ((40]l 
is invariant under the substitutions a —>■ ua and y ^ uy with to = exp(27ri/3). As 
the second transformation leaves the dephased Hadamard matrix invariant, this estab- 
lishes an equivalence between the Hadamard matrices associated with points in D and 
in D' (which one obtains from B through a rotation by 2n/3). As a result, the region ID 
is found to consist of six equivalent sectors, and one may restrict a by 



< arg(a) < 



(42) 
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The transposed SzolWsi family, X^{a, h) is obtained by transposing X(a, h) or by using 
the equivalence H{x,y,u,v)'^ H{x,y,v,u). Fig. [T] illustrates that the points on the 
boundary of the reduced fundamental region for both X{a, h) and X^(a, h) correspond 
to the members of the Hermitean family. 

B Simplification of the Fourier equations in dimension 6 

The conditions for a state \v) G to be MU with respect to Fq are given by P = 
where V = {p± ,q±,r±} with 

p± = —5 ± 2 X5 + 2 X4 ± 2 ^3 + 2 a:2 ± 2 xi + xs^ ± 2 X4X5 + Xi' + 2 X3X5 it 2 X3X4 
+X32 ± 2 X2X5 + 2 X2X4 ± 2 X2X3 + X2^ + 2 X1X5 ± 2 X1X4 + 2 X1X3 ± 2 X1X2 

+2;!^ + 1/5^ ± 2 ^42/5 + ^4^ + 2 7/32/5 ± 2 7/3^4 + ^3^ ± 2 7/22/5 + 2 7/2^4 ± 2 7/22/3 
+2/2^ + 2 2/12/5 ± 2 7/17/4 + 2 7/17/3 ± 2 yi7/2 + 2/1^ , 

q± = -5 + X5 - X4 - 2 X3 - X2 + xi =F v^2/5 =F V^Ui ± >/32/2 ± >/32/i + 2^5^ + a;4a;5 
+Xi' - X3X5 + X3X4 + xs^ - 2 X2X5 - X2X4 + X2X3 + X2^ - X1X5 - 2 X1X4 
-X1X3 + X1X2 + xi^ lb \/3y5X4 ± \/3y5X3 =F s/^y^xi + 7/5^ =F ^37/4X5 ± \/37/4X3 

±a/37/4X2 + 2/42/5 + 2/4^ =F \/3y3X5 T \/37/3X4 ± \/37/3X2 ± \/3y3Xi - 7/3^5 

+2/32/4 + 2/3^ =F \/3y2a;4 =F \/3y2a;3 ± \/37/2a;i - 2 ^22/5 - 2/22/4 + 2/22/3 + 2/2^ 
±\/3yiX5 =F \/37/iX3 =F \/3yiX2 - 2/12/5 - 2 ViVa - 2/12/3 + 2/12/2 + 2/i^ , 
r± = -5 - X5 - X4 + 2 X3 - X2 - xi =F \/37/5 =F \/37/4 ± \/37/2 T VSt/i + xs^ - X4X5 
+Xi' - X3X5 - X3X4 + xs^ + 2 X2X5 - X2X4 - X2X3 + X2^ - X1X5 + 2 X1X4 
-X1X3 - X1X2 + xi^ ± \/3y5X4 =F \/3y5X3 ± \/37/5Xi + =F \/3y4X5 ± \/3y4X3 
^^yiX2 - ym + 2/4^ ± \/3y3X5 =f \/37/3X4 ± V3y3X2 =F ^^37/3X1 - ysys 
-2/32/4 + 2/3^ ± \/3y2a::4 =F V3y2X3 ± V3y2Xi + 2 y22/5 - 2/22/4 - 2/22/3 + 2/2^ 
^ VSyiXs ± VSyiXs =F \/3yiX2 - 2/12/5 + 2 yiy4 - 2/12/3 - 2/12/2 + 2/i^ ■ (43) 

Upon substituting the normalization condition {v\v) = 1, or 

xi^ + yi^ + X2^ + y2^ + X32 + y3^ + X4^ + ^4^ + + ys^ = 5 , (44) 

one finds 

P++P- = , 
p+ - p_ - g+ + g_ + r+ - r_ = , 
2p+ - 2 p_ + (/+ - g_ - r+ + r_ = , 

p+ lb p_ ^ r+ — r_ = , (45) 

giving Eqs. ((16)) . 
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